Detectors for probing relativistic quantum physics beyond perturbation theory 
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We develop a general formalism for a non-perturbative treatment of harmonic particle detectors 
in relativistic quantum field theory using continuous-variables techniques. By means of this we 
forgo perturbation theory altogether and reduce the complete dynamics to a readily solvable set 
of first-order, linear differential equations. The formalism applies unchanged to a wide variety 
of physical setups, including arbitrary detector trajectories, any number of detectors, arbitrary 
time-dependent quadratic couplings, arbitrary Gaussian initial states, and a variety of background 
spacetimes. As a first set of concrete results we prove non-perturbatively, and without invoking 
Bogoliubov transformations, that an accelerated detector in a cavity evolves to a state that is 
very nearly thermal with a temperature proportional to its acceleration, allowing us to discuss the 
universality of the Unruh effect. Additionally we quantitatively analyze the problems of considering 
single-mode approximations in cavity field theory and show the emergence of causal behaviour when 
we include a sufficiently large number of field modes in the analysis. Finally, we analyze how the 
harmonic particle detector can harvest entanglement from the vacuum. We also study the effect 
of noise in time dependent problems introduced by suddenly switching on the interaction versus 
ramping it up slowly (adiabatic activation). 



PACS numbers: 04.62. +v 

I. INTRODUCTION 

For many years, the well-known Unruh-Dewitt model 
^ has been used to explore aspects of quantum field 
theory in curved spacetimes. The great success of this 
model, which couples a qubit to a quantum field using 
a simple monopole interaction, has been its use in an- 
alyzing the observer dependence of relativistic quantum 
phenomena. For example, it has provided satisfactory re- 
sults in the study of phenomena like the Unruh effect [2] , 
meaning that the response of an accelerated qubit detec- 
tor is thermal with the characteristic Unruh temperature. 
This result does not require the use of Bogoliubov trans- 
formations between inequivalent field expansions and the 
subsequent tracing over degrees of freedom beyond a hori- 
zon, but is instead a consequence of a direct calculation 
of the response of the detector when traversing a timelike 
hyperbolic trajectory in spacetime yjj. Additionally, the 
Unruh-Dewitt model is actually a very good basic de- 
scription of the light-matter interaction and reproduces 
quite well the interaction between atoms and light when 
no exchange of angular momentum is involved 4J. 

The main shortcoming of this model is that it is lim- 
ited to perturbation theory. One is therefore barred from 
using it to study problems in which a perturbative expan- 
sion is not a good approximation. These include strong 
coupling, long times and high average energy exchange 
processes. 

We propose here to model a detector as a quantum 
harmonic oscillator rather than a qubit, an idea that has 
been proposed before in other contexts [SHIO]. In other 
words we simply replace two energy levels with infinitely 



many evenly-spaced levels. Nevertheless, qubits are, in 
many cases, just approximations to systems with many 
more levels, so in some ways our description for a parti- 
cle detector is more natural. Given that most symmetric 
potentials in nature can be approximated by a harmonic 
potential for low energies, a harmonic-oscillator detector 
can model a wide range of detectors, from atomic elec- 
tromagnetic levels to the molecular vibrational spectrum. 
In particular, we will consider such detectors in the con- 
text of cavity fields (i.e. the fields they interact with will 
present an IR cutoff), meaning that the field modes are 
discrete. 

Using an oscillator detector has significant advantages 
over the standard Unruh-DeWitt (qubit-based) detector. 
First, the quantum evolution can be solved nonperturba- 
tively. This results from using the symplectic formalism 
for Gaussian states and operations [11 . Many of the sce- 
narios of interest in relativistic quantum theory involve 
quadratic Hamiltonians, making this formalism widely 
applicable. 

Second, the evolution can be evaluated by simply solv- 
ing (in general numerically) a set of coupled, ordinary, 
first-order, linear differential equations. Furthermore, 
the form of this ODE is universal, meaning that one can 
solve a large range of problems with a very minimal ef- 
fort. In particular, this approach can be used to solve 
(a) arbitrary time-dependent trajectories, (b) arbitrary 
quadratic, time-dependent interaction Hamiltonians and 
boundary conditions, (c) arbitrary Gaussian initial states 
of the field modes and detectors, (d) any number of cavity 
modes, and (e) any number of detectors. As we will see, 
a wide range of different scenarios can therefore be solved 
non-perturbatively by the same simple differential equa- 
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tion, which imphes considerable explanatory power and 
computational gain. For example, as we will see, there is 
no need to repeat the numeric calculation whenever we 
want to change a given initial state if the time-dependent 
Hamiltonian mediating the interaction is the same. This 
is rather unlike the perturbative Unruh-DeWitt model in 
which considerably more effort is required. This univer- 
sality, plus the ability to sidestep perturbation theory, is 
the true power of this approach to detector models. We 
will demonstrate this in Section Hvl 

One obvious limitation of this approach is that to solve 
the equations in practice, one is forced to apply an in- 
frared cutoff to the field. However, an infrared cutoff 
naturally appears when studying quantum field theories 
in finite volumes (e.g., optical cavities, periodic waveg- 
uides, etc.), and so this formalism enables us to non- 
perturbatively solve problems of quantum field theories 
in curved spacetimes inside cavities, a matter of great 
interest that has not been thoroughly explored to date. 
If a tabletop experiment in which relativistic quantum 
phenomena is to appear, discrete systems P"2 or super- 
conducting circuits [f3, 14 have an edge on experimental 
feasibility. 

Although in practice one is also forced to use a UV 
cutoff (namely, computing with only a finite number of 
modes), in the results presented in this paper we have 
been careful to find a convergent solution with respect 
to the number of field modes. Specifically, we run the 
simulation with more and more modes until the results 
do not change anymore. As such, this is not a practi- 
cal limitation. These aspects are especially important in 
section [IV B[ where we study explicitly the effects of an 
UV cutoff on the causal structure of our setting. 

The idea of using harmonic oscillators in relativistic 
quantum field theory as particle detectors to obtain non- 
perturbative results was explored by Bei-Lok Hu and col- 
laborators, who reported interesting analytical results in 
[8j. Along with its considerable technical accomplish- 
ments, this approach emphasized that the Unruh effect 
is not reliant on gravitational or geometrical arguments, 
but can be understood as a dynamical effect insofar as it 
indicates how the quantum vacuum affects the response 
of a detector contingent on its motion. In general, a de- 
tector detects field quanta with a nonthermal spectrum, 
where the degree of nonthermality is governed by the pa- 
rameter that measures the deviation from uniform accel- 
eration jl5j . However the practical scope of this approach 
remains to be seen — thus far it has been limited to very 
concrete problems in relativistic quantum theory jl6H18j 
due to their complexity and the number of assumptions 
and approximations required to obtain quantitative re- 
sults. 

In performing our analysis we shall employ a more 
powerful Gaussian formalism, which provides a more ef- 
ficient way to address problems of time evolution when 
considering quadratic Hamiltonian and Gaussian states. 
In this sense our approach is similar to that of Dragan 
and Fuentes, [T^ who made use of the Gaussian formal- 



ism to study a multimode time independent quadratic 
Hamiltonian of two coupled harmonic oscillators. This 
approach had some advantages insofar as it did not re- 
quire any perturbative approximations. However their 
analysis was limited to 1) a single field mode and 2) time- 
independent Hamiltonian. Under that proviso, only sta- 
tionary scenarios and very simple trajectories of detectors 
can be considered. To study a particular non-inertial sce- 
nario (namely eternal uniform acceleration) they relied 
on the existence of Bogoliubov transformations between 
inertial modes and Rindler modes, rendering thermality 
an a-priori assumption instead of a consequence. Fur- 
thermore, by applying free Bogoliubov transformations 
to a single field mode, they were unable to see border ef- 
fects when analyzing the Unruh effect in cavities. Indeed 
the applicability /validity of continuum Bogolibov trans- 
formations for eternally accelerating observers in cavity 
settings in any regime is a rather obscure topic that has 
not been thoroughly understood to date. 

In what follows we present a way to work with an arbi- 
trary time dependent quadratic Hamiltonian and an arbi- 
trary number of modes, being able to analyze scenarios in 
cavities for arbitrary trajectories of an arbitrary number 
of detectors coupled to the field in the cavity without the 
need to assume any Bogoliubov transformations. Fur- 
thermore we can overcome causality violation problems 
[19] of single mode detectors undergoing general trajec- 
tories. We will see that our results are devoid of faster- 
than-light signalling, unlike previous results limited to 
single mode approaches. 

Far from being a mere presentation of some mathe- 
matical tools, we here obtain non-perturbative answers 
in very interesting and yet unstudied cavity scenarios: 

(1) Effects of sudden switching: It is known that the 
response of a detector that is very carefully switched is 
very different compared to the interaction being suddenly 
turned on [2^ . We present a brief non-perturbative anal- 
ysis of this problem showing in a very direct way that it 
is possible to smoothly switch on the interaction without 
exciting the detector. 

(2) Causal signaling inside cavities: It is well known 
that imposing a cutoff in the number of field modes al- 
lows for acausal signaling |19| . This is not surprising 
since a propagating perturbation in a cavity cannot be 
expanded in terms of a finite number of stationary waves. 
By starting one of the detectors in an excited state and 
seeing how long it takes the other one to notice its pres- 
ence we will demonstrate how causality is recovered in 
our setting as we increase the number of modes of the 
cavity. We find the minimum number of cavity modes 
that must be modelled in order to ensure causality for a 
given setup. Note that to completely analyze this kind 
of process perturbatively would require an analysis up to 
fourth order. 

(3) The Unruh effect in a cavity: The question as 
to what the response of an accelerated particle detector 
would be inside a cavity has not been properly explored 
yet. We will show that the response of an accelerated 



detector inside a cavity is still thermal, with some cor- 
rections coming from the border effects of considering 
cavities. Our work provides evidence indicating that the 
Unruh effect occurs not only for the standard Unruh- 
DeWitt model, but also for a discretized harmonic os- 
cillator model. Note that we will not rely on any Bo- 
goliubov transformations or quantization process in any 
accelerated frame. We just answer the question of what 
the response of an accelerated detector is inside a cavity. 
We thereby obtain the Unruh effect inside cavities with 
a derivation similar to the well-known Unruh-DeWitt re- 
sult for the continuum of modes [3] that additionally fully 
sees the appearance of border effects. Furthermore, we 
see that the Unruh effect is also present in a completely 
non-perturbative calculation. The fact that we can eas- 
ily modify trajectories and interaction types means that 
we can easily test the model dependence/independence 
of phenomena like the Unruh effect and entanglement 
harvesting. Such knowledge is very important for under- 
standing the true physicality of such effects in further 
research. 

(4) Vacuum entanglement harvesting: We use the har- 
monic oscillator detector model to analyze a scenario pre- 
viously studied in the standard Unruh-Dewitt setup: the 
extraction of vacuum entanglement by two inertial de- 
tectors [m [22] . We will also comment on the differences 
from the harmonic oscillator model and the qubit model 
in order to see spacelike entanglement. 

II. THE MODEL 

A. Physical setup 

One might suspect that replacing a qubit (two energy 
levels) with an oscillator (infinite energy levels) in our de- 
tector model would complicate the problem, but instead 
it simplifies it, a fact that has been pointed out previ- 
ously [To] . The essential feature that makes this possible 
is that all states in the problem are Gaussian with no 
displacement, and all evolutions are homogeneous Gaus- 
sian unitaries. This means that all states have Wigner 
functions that are Gaussian with zero mean, and all evo- 
lutions are generated by Hamiltonians that are quadratic 
in ladder operators with no linear terms. Such Hamil- 
tonians preserve the Gaussian nature of the states [TT] . 
This means that we don't have to keep track of the en- 
tire Wigner function; all we need to evolve is the co- 
variance matrix, whose size scales quadratically with the 
number of modes (rather than exponentially, as is the 
case for general states). Similarly, all evolutions can be 
represented by symplectic matrices, which has the same 
scahng [53]. 

For calculational purposes, we assume that an IR cut- 
off of some length L has been imposed on the fieldj^ As 



This is necessary because we want to use matrix algebra to nu- 



3 




FIG. 1. Harmonic-oscillator detector (red dot) moving 
through a fixed cavity. This is to illustrate the difference be- 
tween our setup and those in which the cavity itself is in mo- 
tion [24]. Note that any number of detectors may be present. 
This particular case of an accelerated detector is treated in 
Sec UVCl 

such, our physical model is that of a detector moving 
around in a large cavity. There is an important distinc- 
tion to be made here with other models that consider 
the cavity itself as being in motion [24'. In our work, by 
contrast, the cavity is large and fixed, and the detector 
moves within it (if it moves at ah). See Fig.[l] 

The Hamiltonian that we propose below is very sim- 
ilar to one already discussed in the literature in which 
a (not fully general) Hamiltonian that describes the in- 
teraction of a harmonic-oscillator detector with a finite 
number of field modes was presented [10] . Although they 
focus was on an approximation of time-independent cou- 
pling to just a single field mode (a special case of our ap- 
proach), the advantages of the Gaussian formalism and 
the generality of this approach was made evident. The 
time-independent single-mode approximation |10| has the 
advantage of simplifying the obtention of analytic solu- 
tions but at a price of relying on Bogoliubov transforma- 
tions for uniformly accelerated detectors and limiting to 
regimes where a single mode approximation can be valid, 
whereas our results are generally numeric but arise from 
an approach that has much greater applicability. 

B. Hamiltonians generating evolution with respect 
to different time parameters 

In relativistic scenarios it is important to keep in mind 
that Hamiltonians generate evolution with respect to a 
given time parameter that does not necessarily coincide 
with the proper time of some (or any) of the proper times 
of the physical subsystems that are in interaction. When 
we consider the global Hamiltonian of multipartite sys- 
tems we would need to express it in terms of a common 
time parametrization. Due to this, we also wish to pro- 
vide a discussion of how to generate evolution with re- 
spect to an arbitrary time parameter. In particular it 



merically solve the resulting differential equations, altfiough a 
formal generalization of our metfiod to the continuum limit may 
be possible in the form of integro-differential equations resem- 
bling those in Section |III[ A complete formulation of this is left 
to future work. 
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can be in general useful to evolve in a global time coordi- 
nate t, particularly in the case where there are multiple 
detectors in which each detector j has a different proper 
time coordinate Tj associated with it. The calculation is 
most straightforward in the Heisenberg picture, although 
applies equally well in the Schrodinger or interaction pic- 
tures. In this way we provide a "dictionary" by which we 
can transform to any other time coordinate. In Sec. |II C| 
we continue by introducing the general form of Hamilto- 
nians that can be used with our approach. 

Let us proceed, then, by considering a general time- 
dependent Hamiltonian H(t), which generates transla- 
tions of the entire system in the global time coordi- 
nate t. This Hamiltonian includes the free Hamiltonian 
for each system, as well as interactions. With respect 
to t, the Heisenberg equation of motion for a general op- 
erator A{t), possibly having explicit dependence on t, is 



' Ait) ^ Urn, Ait)] + 



dt 



dt 



(1) 



A different choice of time coordinate can be taken into 
account by applying the chain rule. For the moment, let 
us choose the new time variable to be the local proper 
time Tj that parametrizes the world line [xirj), ^(tj)] tra- 
versed by detector j. Applying the chain rule gives 



d I r / NT dt d ~ , . 



t = t{Tj) 

dt dAit) 



dt 

dTj 



H[tiT,)]),A[tir,)] 



, dA[tiT,)] 



(2) 



Thus, we can start with the Hamiltonian Hit), which 
generates translations in the global time coordinate i, 
and then define 



4.(r,) := ^H[tir,)] 



(3) 



as the Hamiltonian (for the entire system) as seen by 
detector j, which generates evolution for the entire sys- 
tem with respect to the proper time coordinate Tj. The 
derivative dt/drj is the redshift factor for an observer in 
the detector's reference frame, which provides an overall 
scaling of all energies in the combined system (because 
this is what such an observer would experience). Notice 
that although the notion of proper time is local, we need 
to be able to evolve the entire system with respect to 
this coordinate because we are working in the Heisen- 
berg picture. This is not a problem as long as ^(tj) is 
an invertible function over the range of times of interest. 
The equivalence of the two pictures is made explicit by 
defining Heisenberg operators that arc more natural to 
the detector's frame: 



Mr,) := A[tirj)] 



(4) 



We can now use Eqs. ^ and Q to rewrite Eq. ^ as 



dr. 



(5) 



For example, if Ait) were to represent the position of 
the second hand on a wristwatch worn by an observer 
traveling with detector j, then it would make more sense 
to consider Ajirj) because this operator would have a 
simpler evolution with respect to Tj than Ait) would with 
respect to t (since the wristwatch evolves more simply 
with respect to Tj than with respect to t). Similarly, 
it will be easier to start with the simple version of the 
wristwatch's Hamiltonian HjiTj) and then invert Eq. ([s]) 
to obtain 



Hit) 



dT. 



dt ^' 



(6) 



(which will be more complicated) for use in the global 
Hamiltonian. 

The upshot of all of this is that we can define a single 
Hamiltonian Hit) for the whole system with respect to 
some global time coordinate t and then use Eq. ^ to 
transform it to any other time coordinate we wish to 
use for the evolution. Furthermore, when building up 
this Hamiltonian, it will sometimes be easier to start by 
defining a piece of it with respect to local proper time 
and then use Eq. Q to figure out what this piece looks 
like in the global time coordinate. 



C. The Unruh-DeWitt Hamiltonian in general 
scenarios 

We have to be careful when we want to deal with 
Hamiltonians generating translations with respect to dif- 
ferent time parameters above all when we want to de- 
scribe the interaction of systems that have different 
proper times. 

Indeed, in general scenarios it is not trivial to define 
either the interaction nor the free Hamiltonian in dif- 
ferent pictures. To guide the reader through this sec- 
tion let us introduce the following notation: We will call 
H^ , H^ , H^ respectively the complete Hamiltonian in 
the Schrodinger, interaction (Dirac) and Heisenberg pic- 
tures. We will denote with the subscript '0' the free part 
of the Hamiltonian and '1' the interaction part. 

Also we will include a superindex t or t denoting with 
respect to which time the Hamiltonian is a generator of 
translations. 

We will consider the interaction of a number of particle 
detectors with a quantum field. To model this interac- 
tion we will consider an X-X coupling of the form of the 
Unruh-DeWitt Hamiltonian [T]. Note however that the 
formalism we present is much more general than this and 
we can in fact use any quadratic Hamiltonian that we 
like. For out immediate purposes we choose to use the 
X-X coupling in order to compare with previous works. 



5 



For every detector coupled to the field the Unruh- 
DeWitt interaction will be described by the following 
Hamiltonian 



i/i =A(T)A0[x(r)], 



(7) 



where A(r) is the switching function, fi is the monopole 
moment of the detector and (f>[x{T)] is the field op- 
erator evaluated along the worldline of the detector 
parametrized in terms of the time r with respect to which 
the Hamiltonian generates translations . 

For the sake of clarity let us start our reasoning with 
a very simple scenario: let us consider a single detector 
undergoing general motion in flat spacetime with an as- 
sociated proper time t and a scalar quantum field that 
we will choose to expand in terms of plane wave solu- 
tions in terms of a global Minkowskian time t, as it is 
commonplace in quantum field theory. 

To derive the correct form of the Hamiltonian in the 
interaction and Heisenberg pictures let us write first the 
field and monopole operators in the Schrodingcr picture: 



n 

fL^ = {ad + flrf), 



(8) 
(9) 



where Vn{x) are the spatial part of the solutions to the 
field equations, i.e. the mode functions. For instance in 
the case of reflecting boundary conditions these would 
be Vn{x) = sin(A:„a;) whereas in the case of periodic ones 



Vnix) 



. Here k„ is the wavevector for field mode 



n; we will specify its form below. 

Let us start from a very well known result from first 
principles: we can write the the free Hamiltonian for the 
field, and the free Hamiltonian of the detector in their 
respective times in the Schrodingcr picture 



frS,t 

-"O, field 



E 



LO, 



1 ttj, 



h: 



O.dot 



(10) 
(11) 



Now, to write the complete free Hamiltonian we can- 
not just naively sum these two terms together because 
they generate translations with respect to different time 
parameters. We would need first to transform them to a 
common time parameterization. We will see that in or- 
der to recover the correct form of the well-known Unruh- 
Dewitt Hamiltonian in the interaction picture, we must 
transform the field Hamiltonian to generate translations 
in the proper time of the detectors. In this way, using 
^ we have that 



0, field 



d 

d7 



(12) 



so that we can write the complete Hamiltonian in the 
Schrodingcr picture generating translations in r as 



h; 



S.T 



H 



S,r 



being 



(13) 



n 

(14) 

Note that the free Hamiltonian is not time independent 
as it was in the case where the detector is inertial. 

In most textbooks [3], calculations involving non- 
inertial detectors coupled to the field are dealt with in 
the interaction picture. We will see that we recover the 
well-known form of the interaction Unruh-DeWitt Hamil- 
tonian by changing from the Schrodingcr to the interac- 
tion picture. 

Recall the transformation between the Schrodingcr and 
the Interaction pictures: 



H 



D,T 



C/T(r)ff^^^C/o(r), 



(15) 



where Uq{t) is the solution to the Schrodingcr equation 



in r using just the free Hamiltonian Hq'^ : 



(16) 



Notice that in this case the transformation is non-trivial 
due to the non-trivial dependence on r of the global time 
parameter t(r). This yields an explicit time dependence 
of the field free Hamiltonian. Since Hn'^ commutes with 



itself at different times, we can solve Eq. ( 16 1 explicitly 



without needing to worry about time ordering: 

Uo{t) — exp —i dr h 
V Jq 

exp ■ ' 
exp 



n 

i( ^'^naJian)i(T) - if^ajjadT , (17) 



This operator leaves invariant the free parts of the Hamil- 
tonian, and its action on the a„ and operators is 



Ul{T)aMT) = e 



_ -iw„t(r)- 



Ul{T)ddUo{T) ^ e-'''^ ad 



(18) 
(19) 



allowing us to write the Unruh-Dewitt Hamiltonian in 
the interaction picture with respect to the parameter t 
as 



H 



: a;„4a„ -f na\dd + {dde '^^^ -|- d\e'^^) 

n 

Kt) E [anUn[x{T), t(r)] -I- 4<[x(t), i(r)] , 

n 

(20) 
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where 



(21) 



Thus we recover the standard form Unruh-Dewitt Hamil- 
tonian [3] in the interaction picture that generates trans- 
lations with respect to time r starting from the well 
known free Hamiltonians (in the Schrodinger pictures, 
with respect to their respective natural time parameters) 
after transforming to a common time r and changing to 
the interaction picture. 

Notice that the (real-valued) coupling parameters A(t) 
can be time dependent. This allows for, among other 
things, switching the detector on and off with a particu- 
lar temporal profile known as the switching function or 
switching function. The last pair of terms indicates that 
the detector interacts with the field with a strength (and 
phase) governed by the mode functions m„(x, t). In order 
to determine these mode functions, we need to establish 
boundary conditions for the cavity. These can include, 
for example, Dirichlet boundary conditions, in which the 
field strength vanishes at the boundary (as in a physi- 
cally realistic optical cavity) or periodic conditions (as 
in a physically realistic periodic waveguide). The mode 
functions for these cases take the forms 



/ .N Dirichlet / . ,\ ■ n \ 

Un\X, t) I > exp ( — IWntj sm {knX) 



and 



u„{x,t) 



% exp (— iw„i) exp (iknx) , 



(22) 



(23) 



where, as we are going to work with massless fields, w„ = 
|fc„| and kn = nn/L or fc„ = 2mr/L respectively for the 
Dirichlet or periodic boundary conditions. 

Now, we will find it convenient for use in the next 
section to work in the Heisenberg picture. We will use 
the fact that the form of the complete Hamiltonian in 
the Heisenberg picture coincides with the form of the 
Hamiltonian in the Schrodinger picture. To see this, note 
that the transformation between the two pictures is by 
the full time-evolution operator C/(r), which satisfies the 
full Schrodinger equation 



dr 



(24) 



The Hamiltonian in the Heisenberg picture is obtained 
from its Schrodinger-picture counterpart by the usual 
transformation between the two pictures for any oper- 
ator: 



H 



H,r 



(25) 



This means that we do not have to do any work to mod- 
ify the Schrodinger-picture Hamiltonian in order to use 
it in the Heisenberg picture. All we have to do is reinter- 
pret all operators within it as being Heisenberg-picture 
operators instead of Schrodinger-picture ones. 



After all these simple steps, it is very easy to write 
the most general X-X type Hamiltonian for an arbi- 
trary number of detectors undergoing general trajectories 
with different proper times Tj and with time dependent 
couplings. However, if multiple detectors have different 
proper times then we again need to be careful. One must 
always make a choice of time, and in this more general 
case it makes sense to use the global Minkowski time t. 
Transforming the Hamiltonian to time t and reinterpret- 
ing all operators in the Heisenberg picture yields 

N M J , s 



ji=i 

N 



+ + 4 J [anVnixjit)] + aivn[xj{t)] 

n=l 

(26) 

where Xj{t) is the trajectory of the j-th detector 
parametrized in terms of the global Minkowskian time 
t, and all operators are now understood to be in their 
Heisenberg representation. We will see in the next section 
how working in this representation allows us to derive a 
simple, number-valued equation of motion that describes 
the full evolution of the detectors-|-field state. 

To recapitulate, the Hamiltonian above represents a set 
of -t- M time-dependent coupled harmonic oscillators. 
M of them are oscillator-based Unruh-DeWitt detector 
modes (labeled with the index dj), and the other N are 
modes of the quantum field inside a cavity (labeled with 
an integer index n). Notice that the there is no direct 
detector-detector coupling, and the field is a free field, 
meaning there is no coupling directly between field modes 
either. 



We emphasize that Eq. (26) is not the most general 



form of the Hamiltonian that could be imagined in this 
scenario. It is simply the same as the original Unruh- 
DeWitt detector model. The connection is made by 
choosing no relative phase between ddj and aj^ in the 
interaction term, which makes their sum proportional to 
the position (monopole moment) of the oscillator. In 
general our formalism, to be presented now, is trivially 
capable of solving far more general interaction models, 
provided they are quadratic. 

Our problem is now this: given detector worldlines 
[t{Tj) , x{Tj)] and an initial (Gaussian) state for the detec- 
tors and field, evolve the detectors and the field using the 
Hamiltonian in Eq. ( 26 ) , and consider the reduced state 
of the detectors after the evolution. In order to make 
use of the simplification afforded by the use of Gaussian 
states and quadratic Hamiltonians, in the next section 
we will derive a differential equation using the symplectic 
formalism for the Heisenberg picture and another differ- 
ential equation using the Hilbert space evolution in the 
Interaction picture. Working with either of these pictures 
will let us compute the covariance matrix for the field and 
detectors throughout the evolution, and since the state 
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remains Gaussian the whole time, this is equivalent to 
tracking the evolution of the full state itself. 



III. SOLVING THE EVOLUTION 

In the case of an oscillator detector the great advan- 
tage is that one is able to utilize the Gaussian formalism 
[m [231 HSl [26]. In this solution method we will find it 
convenient to use the Heisenberg picture. As discussed in 
Section p^I C[ there are subtleties involved because of the 
different time coordinates available (proper time for each 
detector, plus the global time coordinate of the field). 

In section [III Aj we will use the Heisenberg-picture 
Hamiltonian, which generates translations in the global 
coordinate time t, to derive a simple linear differential 
equation using the symplectic formalism for Gaussian 
evolutions. This is the simplest, most straightforward, 
and computationally most efhcient approach. 

In section jlll B[ we describe another method that di- 
rectly calculates using the interaction Hamiltonian in the 
interaction picture. This method is more complicated 
than the symplectic method and yields a non-linear sys- 
tem of ODEs, but because it is independently derived, 
and because the resulting differential equations are com- 
pletely different (second-order instead of first-order), this 
provides an independent check of our numerical results. 

The purpose of presenting in parallel both approaches 
is dual: on the one hand we can tackle problems with 
two different approaches (which provides a way to opti- 
mize our computational strategy in order to solve a given 
problem). On the other hand we see that the same results 
can be obtained via independent methods. This serves 
as a connection between the two formalisms. 



A. Phase space evolution 



We will use the Hamiltonian (26) in the Heisen- 



berg picture to evolve Heisenberg quadrature opera- 
tors {qdj{t),pdj{t)) for each detector and {qn{t) , Pn{t)) 
for each field mode. To streamline calculations, we stack 
these operators on top of each other to form the fol- 
lowing vector of operators, while omitting the explicit 
t-dependence for clarity: 



X := ((7di, ■ 
where 



,qN,Pdi,---,PdM^Pi^ 



,Pn) 



V2 



(ai + al), = — (a^-a) 



(27) 



(28) 



are respectively the canonical position and momentum 
of every single oscillator in the Heisenberg picture. Note 
that the transpose operation merely transposes the 
shape of an operator-valued vector; it does nothing to 
the operators themselves. 



A Gaussian state is one that can be expressed purely in 
terms of its first and second quadrature moments. If we 
neglect phase-space displacements, this means that such 
a state is fully characterized by a covariance matrix cr, 
the entries of which are 



(29) 



Specifically, the Wigner function for the zero-mean Gaus- 
sian state associated with cr is 

I^(x) = 7r-(^+*^)(dct(T)-i/2exp(-xT,T-ix), (30) 

where x is the vector of c-number- valued coordinates cor- 
responding to x: 



{qdi,- ■ ■ , qdM,qi, ■ ■ ■ ,qN,Pdi, ■ ■ ■ ,PdM^Pi^ ■ 



,pn) 



(31) 

Although in general, displacements will be required for 
a full description of Gaussian states and evolution, if 
we start with a zero-mean initial field state (such as 
the Minkowski vacuum) and a zero-mean initial detec- 
tor state as well, then the lack of linear terms in our 
Hamiltonian means that the state remains at zero mean 
at all times. 

To compute the evolution, we utilize the fact that 
quadratic Hamiltonians preserve Gaussianity jllj . This 
means that quadrature operators get mapped to linear 
combinations of quadrature operators: 



t' = U^iU = Si, 



(32) 



where C/ is a Gaussian unitary (i.e., a unitary transfor- 
mation generated by a quadratic Hamiltonian). In this 
equation, S is a symplectic matrix of c-numbers that acts 
via matrix multiplication on x as a vector, while U is a, 
unitary operator that acts on the individual operators 
within X. Specifically, 



U^XjU 



2(N+M) 
k=l 



(33) 



Notice that in general there would be a phase-space dis- 
placement term, which would give x' = Sx -|- y, but we 
are neglecting this as justified above. The symplectic 
nature of S is guaranteed because the commutation re- 
lations must be preserved, giving rise to a symplectic 
form n to be preserved by the Heisenberg matrix ac- 
tion. The explicit form of may be deduced by writing 
out the commutation relations for x and requiring them 
to be unchanged under the Gaussian unitary operation. 
Using the notation of Ref. [5^] , the canonical commuta- 
tion relations [qj,pk] = iSjk (with H = 1) can be written 
succinctly as 



I 
I 



=: in, 



(34) 



where the commutator of two operator-valued vectors is 
defined as 



[f,s'^] := fs'^ - (sf'^)T . 



(35) 



8 



Therefore, ft has elements fli 



-i[xi,Xj'\. That S is 



symplectic means that SJIS"^ = ft, which follows from 
requiring that the new commutation relations must equal 
the old ones: [Sic, (Sx)"^] = iCl. (See Ref. ^B] for more 
details.) 

A general quadratic Hamiltonian generates a Gaussian 
unitary U(t) = g"*/^"'*''* that is associated, by Eq. ( [32| , 
with a symplectic matrix S{t), which satisfies 



x(i) = C/(t)TxoC/(t) = S(i)xo , 



(36) 



where all time dependence (or lack thereof) is indicated 
explicitly, and xg is the initial vector of quadratures at 
t — 0. Correspondingly, the Schrodinger evolution of the 
state, as given by the evolution of the covariancc matrix, 
takes the form 



rrit) ^ S(t)rToS(t)^ 



(37) 



where ctq is the initial state. Our goal in this section is 
thus to find a differential equation for S{t) that represents 
the evolution generated by a quadratic Hamiltonian. 

A general time-dependent, quadratic, Heisenberg- 
picture Hamiltonian H with generated time coordinate t 
can be written as 



H = x'^F(t)x, 



(38) 



where F(t) is a Hermitian matrix of c- numbers containing 
any explicit time-dependence of the Hamiltonian. We can 
now write the Heisenberg equation for the time evolution 
of the quadratures: 



dt 



X — i [H, x] 



(39) 



Writing this out in components and using Eqs. (34) 
and (38 1 gives 



df 



mn 

(t) ^X^n^jn ^jm^n 



which can be collected back into vector form as 
d 



where F*^ 



dt 

(F + F^ 



x = nF''>""(t)x, 



(40) 



(41) 



We now plug in Eq. ( 36 ) , giving 

|[S(t)]xo = f2F^^'°(t)S(t)xo, (42) 

where we used the fact that xq is time independent. Now 
we use Eqs. (34) and (35) to eliminate the operators by 



taking commutators with Xq on both sides: 
d 



dt 



S{t) ] xo,xJ 



|s(t)[xo,iJ] 



fiF^^-(t)S(<)Ao,xJ 



f2F^>'"Xi)S(i)[io,io] 



(43) 



where we have factored out the C-number matrices multi- 
plying Xq [26]. Since [xg, Xq ] — in, which is an invertible 
matrix of c-numbers, we can cancel it, yielding the fol- 
lowing first-order, linear, ordinary differential equation 
for the symplectic matrix: 



dt 



S{t) = f2F"y'"(t)S(t). 



(44) 



Solving this equation with the initial condition S(0) = 
I such that Xq = S(0)xo is completely equivalent to 
solving the standard Hilbert space evolution with the 
Hamiltonian-unitary formalism after taking advantage of 
the quadratic nature of the Hamiltonian, as we will see 
in the next section. 

It will be convenient to know the form of the Hamilto- 
nian in terms of the annihilation and creation operators 
of the system, since this is how the monopole-monopole 
coupling is typically given. To this end, let us stack lad- 
der operators on top of each other to form the following 
column vectors: 



a (fldi , • • ■ , fldM , Si , ■ • • Oat) , 

:= {d\^,...,dl^^,a\,...d^^f. (45) 

The Hamiltonian from Eq. ( |38[ ) can be put into the form 

H = {k^fw{t)a + {a^fg{t)k^ + aTg(t)"a , (46) 

where w(t) and g{t) are coefficient matrices, and = 
M*'^ — M.^* denotes the Hermitian conjugate of any 
matrix M. 



By equating ( 46 ) to ( 38 ) and comparing coefficients it 



is easy to obtain the form of the matrix F(t), which takes 
the following block form: 



A{t) X(t) 
X(t)H B(t) 



where 



Ait) 
B{t) 
X{t) 



(w(t) 
(w(t) 



2 (w(t) 



lit) 
lit) 
lit) 



(47) 

(48) 
(49) 
(50) 



Note that in the simple particular cases where the 
Hamiltonian at different times commute the solution 
of ( [44| ) can be obtained analytically and it is simply 
S(i) = exp(nF"y"^t). 

As a final note, we would like to point out that we de- 
rived these results using global Minkowski time t because 
it is the least biased time coordinate when there are mul- 
tiple detectors involved. But any time coordinate can be 
used instead — such as some particular detector's proper 
time — as long as the appropriate transformation is made 
to the Hamiltonian via Eq. Doing so would define 
a new differential equation of the same form as Eq. ( 44 ) 



but evolving in the new time coordinate instead of in t. 
In fact, this is the most direct way to calculate detector 
responses, and it is the method that we use in Section IV 
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B. Hilbert space evolution 

Although the formalism presented in the section above 
is elegant and fully general for quadratic Hamiltonians, 
it is useful to compare our results with an independent 
method. This section outlines the supplementary method 
that we used for this purpose. Without invoking the 
full machinery of symplectic transformations, we can still 
take advantage of the quadratic nature of the interac- 
tion to derive the same results by numerically calculat- 
ing the unitary time evolution operator, an approach that 
is more standard within quantum field theory. Further- 
more we can us this method directly in the interaction 
picture, rather than requiring one to first transform to 
the Heisenberg picture. Here we give only an outline of 
the method; details can be found in Appendix \K\ 

For a quadratic Hamiltonian, time evolution can be 
expressed in terms of displacements, squeezing, and ro- 
tation unitary operations |27j . In particular, the unitary 
evolution we are looking to solve for can then be put into 
the form 

u{t) = ^^^*^s{z{t))Dm))R{m) , (51) 

where 7 is number valued and the squeezing, rotation 
and displacement operators are respectively defined as 



Dif3) = 



(52) 
(53) 
(54) 



where z and (p are matrices, and /3 is a column vector. 
The notation is used to represent the Hermitian con- 
jugate of z, the elements of which are number valued. 
Note that we without loss of generality we can consider 
z to be symmetric because it is only the symmetric part 
that contributes to S{z). Also note that 4> niust be a 
Hermitian matrix to ensure unitarity of R{4>). 

The exact form of these transformations can be ob- 
tained nonperturbatively by employing a technique in- 
troduced by Heffner and Louisell [23. For the cases of 



interests (for the family of Hamiltonians (20)), and as it 
is thoroughly detailed in appendix]^ due to the algebraic 
nature of the interaction, the evolution can be reduced 
to 



U = e'''S{z)R{ct)), 



(55) 



and the problem of solving the dynamics can then be 
reduced to solving the following system of coupled differ- 
ential equations 

iC{t) = 4C,(i)gH(t)C,(t) + 2w(t)C,(t) + g(t), (56) 
it>{t) = (4C,(t)gH(t) + w(t))(D(t) + I), (57) 

where w(t) and g(t) are the Hamiltonian coefhcient ma- 
trices as defined in Eq. ( 46 1 , evaluated in any picture one 



chooses (for example this can be done directly with the 
interaction Hamiltonian in the interaction picture). We 
have defined Cg = (C -I- C^)/2, where the matrices 
and D are identified with the squeezing and rotation by 



C, = itanh(r)e'^ 
D + I = sech(r)e*'^. 



(58) 
(59) 



where z = re'^ = e^^ . Therefore, numerically solv- 
ing the equations (56), (57) with the initial conditions 
C(0) = D(0) = we can non-perturbatively solve the 
time evolution. 

The covariance matrix a of the detector-|-field state is 
defined as in the previous section. Once we have solved 
for the squeezing and rotation operators, in the sense 
that we have solved for z{t) and 4>{t), it is then straight- 
forward to compute the evolved covariance matrix |28j . 
If we split the covariance matrix into the block form 



qq " 
T 

^qp ^: 



(60) 



then straightforward algebraic operations (See Appendix 
|A) Eqs. (|A6]) and (|A9])) allow one to determine the form 
of these blocks. For example if our detector-|-field state 
starts in the vacuum state then the evolved state is simply 
given by a multi-mode squeezed state S'(z) |0), since the 
vacuum is invariance under rotations. In this case the 
blocks of the evolved covariance matrix take the form 



1 



(cosh(2r) + sinh(2r)e'^ + cosh(2r^) (61) 



■sinh(2r'^)e" 



-(cosh(2r) - sinh(2r)e''' + cosh(2r^) (62) 



qp 



-sinh(2rT)e-'^^), 

(cosh(2r) - sinh(2r)e'^ - cosh(2r'^) 



(63) 



-Ksinh(2r'r)e-^^^). 



Once this is obtained then computing the covariance 
matrix for the detector(s) alone is trivial: one must sim- 
ply isolate the rows and columns of cr corresponding to 
the detector modes. Both methods give the same time 
evolution and the same covariance matrix applied to the 
same scenarios, and we used them both independently to 
check our calculations. 



IV. RESULTS 

As discussed in the introduction, we will see below 
how these tools can be used to observe relativistic quan- 
tum phenomena. After examining effects due to sharp 
switching [201 HH ISO] we will go on to study the emer- 
gence of causal signaling with increased number of field 
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modes [TH], the Unruh effect and vacuum entangle- 
ment harvesting [21] |22] . We will be considering these 
effects in the context of optical cavities. Cavities are in 
fact excellent systems to study since experimental physi- 
cists have developed tools for the precise production and 
control of optical states inside cavities, superconducting 
SQUIDs, or microwave guides, and therefore such sys- 
tems are likely to be key in the experimental verification 
of relativistic phenomena [H] . 

In all of the scenarios we consider here we will simplify 
the calculations by using the detector's proper time r as 
our preferred time coordinate, as is done in the major- 
ity of the literature concerning Unruh-DeWitt detectors. 
Doing this is acceptable here because all of the cases we 
consider either involve only a single detector or two de- 
tectors that share the same proper time. In particular 
we will take for our interaction Hamiltonian the usual 
Unruh-DeWitt interaction, namely monopole-monopole 
coupling, but of course with the usual qubit operators 
replaced by their corresponding oscillator operators. 

In the interaction picture the interaction Hamiltonian 
is therefore 

3 

X ^{Un[x{T),t{T)]an + U*Jx{T),t{T)]al), 
n 

(64) 



where it„ is given by either Eq. (22) or (23), depend- 
ing on the boundary conditions we impose. Since it is 
numerically simpler to change to the Heisenberg picture 
and use the equation of motion as presented in section 
|III A[ we use this method for all calculations and use a 
second numerical method based on the formalism pre- 



sented in section III B as an independent check of our 
results. 

As a note, it is important to realize that some of the fol- 
lowing scenarios are actually solved exactly analytically 
(not numerically). We can do this in the cases where we 
have stationary detectors with constant (sharp) switching 
functions since in these cases the Hamiltonian at different 
times commute. Equivalently, in these cases the solution 
to Eq. (pl is simply S(r) = expiflF^'^'^T). 



A. Stimulation of noise due to sudden switching 

Before examining more involved settings such as the 
effect of non-inertial motion on our detector or the vac- 
uum entanglement harvesting, let us discuss first what 
happens for a single inertial detector in a vacuum back- 
ground. As is known, the use of sufRciently sharp switch- 
ing functions A(t) can stimulate excitation of the detec- 
tor even when it is inertial |30j . This stimulated vacuum 
fluctuations can be reduced by increasing the interaction 
time and using smooth switching functions, for exam- 
ple Gaussian time profiles. Indeed we see non-negligible 



inertial excitation when the time of integration is small 
enough, but we will confirm that this vanishes for long, 
Gaussian switching functions as it corresponds to the de- 
tection of quantum noise in the vacuum state of a quan- 
tum field. 

There have been some studies on the effect of the 
smoothness of the switching function on the probability 
of excitation of the Unruh-DeWitt detector and the stim- 
ulation of quantum noise [201 HI] ■ These studies showed 
that the quantum noise that a detector will observe for 
short interaction times is strongly influenced by the way 
in which the detector is switched on. These studies were 
done within perturbation theory and only for the case 
of fields in free space. It would be interesting to show 
these effects in cavity settings and in a non-perturbative 
calculation. 

To this end we consider a single detector sitting in- 
ertially in the center of a cavity of length L, such that 
^(t) = r and x{t) — L/2. Seeing as we are considering an 
optical cavity we will use reflecting boundary conditions 
for the field. We then solve for the evolution generated by 
the Hamiltonian (46) with some switching function A(t). 



After the evolution, our detect or -(-field acquires a multi- 
mode squeezed, pure-state covariance matrix cr. The 2x2 
covariance matrix cr^ corresponding to the oscillator de- 
tector is then obtained by taking the detector-detector 
elements of <r: 



CTd = 



'-' aa ^ a 



'qq 

Ad) 
'qp 



'IP 
Jd) 



(65) 



The symplectic eigenvalue v of this state can easily be 
computed as the absolute value of either of the eigenval- 
ues of the matrix zOcr^ (they come in a ± pair) where $7 
is the symplectic form 



n = 



1 

-1 



(66) 



The value of v gives the mixedness of the state, with 
I' — 1 corresponding to a pure state. To be precise, the 
purity is given by Trp^ = In general it is somewhat 
nontrivial to compute the excitation probabilities p„ = 
{n\^ Pd\n) ^ of a given Gaussian state [SI], but for our 
purposes we will find it sufficient to only consider the 
probability of no excitation, pq. In the case of a zero- 
mean state (i.e. the Gaussian Wigner function is centered 
at the origin of phase space), which includes the states we 
consider (as shown in the appendix [A|, this probability 
is given by [31] 



Po 



dctcr^ + TrcTfl 



1. 



(67) 



To examine the stimulated excitation of the detector 
due to the activation of the interaction we will consider 
two different switching functions. The first will represent 
sharp switching on and off, in which we set the switching 
function to a constant A(r) = A and track the evolution 
through time. In the second case we will use a Gaussian 
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FIG. 2. Excitation probability of an inertial detector: In Fig. (a) we consider a sharp switching function and track the 
excitation probability though time r. Even when r becomes much larger than the values shown in the plot, the probability 
remains nonzero. In Fig. (b) we consider Gaussian switching functions of standard deviation 5 and we plot the excitation 
probability after the detector has finished evolving (i.e., when the Gaussian tail has become negligible) hence 5 in that plot 
is a measure of smoothness. For sharp switching we observe excitation as expected and for Gaussian switching we observe 
excitation only when the Gaussian is sharp enough, which is also as expected. The parameters used are A = 1/00, L = 27r, 
Q, = 9/2 (resonant with the 9th mode) and the detector is placed at position x — L/2 = tt. 



switching function of the form A(r) = Aexp(— t^/2^^) 
and will integrate from time Ti ~ —45 to rj = 45. b is in 
this case a measure of the smoothness of the time pro- 
file. Note that for sharp switching considerably more field 
modes must be included before we observe solution con- 
vergence than in the case of Gaussian switching. This is 
expected since a sharp A(t) can excite field modes signif- 
icantly higher than those near resonance. This is because 
the off-resonant rotating wave terms become important 
(as well as the counter-rotating wave terms) if the in- 
teraction suddenly changes in characteristic times faster 
than ~ 1/il; see [2D]. We then use Eq. (67) to compute 
the probability of excitation, 1 — po: for both the sharp 
switching as a function of r and Gaussian switching as a 
function of 8. The results are plotted in Fig. [2j 

In both cases we see that the excitation probability 
tends to zero as the interaction time goes to zero, as it 
must since this is the limit of no evolution. Note that 
these results can be also easily computed perturbatively, 
giving the same answer (up to higher than second order 
corrections) as that shown in Fig. [2] For the Gaussian 
switching function we see that the excitation becomes 
negligible for larger 5; this results from the switching 
function becoming smoother with increasing b and is ex- 
actly what should be expected. For sharp switching how- 
ever we see that the excitation probability does not decay 
with increasing time because in this case it is the initial 
discontinuity in A(t) that causes the excitation. 



B. Emergence of Causal Signaling 

The causal behaviour of the probability of excitation of 
spacelike separated atoms have been thoroughly studied 
in the continuum in the context of the so-called Fermi 
problem [351 [33] , where satisfactory answers have been 



provided (see among others [131 13i] '). 

Regarding IR cutoff settings, it has been pointed out 
that a single mode approximation in the Unruh-Dewitt 
model (namely, considering a system of detectors inter- 
acting only with a single mode of the quantum field) leads 
to superluminal signaling |19| . This is not surprising: a 
complete set of solutions to the field equations inside a 
cavity are the stationary waves inside it. A propagating 
signal cannot be expressed in terms of a finite number 
of stationary waves. Strictly speaking, to completely re- 
cover causality one should consider the infinite number 
of modes inside the cavity. 

However, it is also known from quantum optics [35] 
that the Jaynes-Cummings model (basically a single- 
mode-approximated Unruh-DeWitt model) produce ac- 
curate results if the evolution times are long. This is so 
because for long times (much longer than the light cross- 
ing time of the cavity) a stationary regime is reached and 
for infinite times there is, of course, no signal propagation 
issues. Hence if we are to study quantum information- 
related topics with this models we need to take seriously 
the issue of causality. This is one reason why, for cer- 
tain scenarios, the use of only a single mode is highly 
unphysical. 

Here we point out the fact that, using our work, we 
can very easily observe the emergence of causal signal- 
ing by considering our detector to be coupled to different 
numbers of field modes. Although rigorously speaking, 
we would need to add up the infinite number of modes 
to completely guarantee causality, we will discuss that, 
for a given arbitrarily high time precision, it is possible 
to find an effective model with a finite number of modes 
such that causal signaling does not happen within the 
required time precision. To do so we will consider two 
inertial detectors in the same cavity, and our goal is to 
examine how long it takes one detector to observe the ef- 
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FIG. 3. Excitation probability as a function of time of an inertial detector (detector 1) in the presence of another in a highly 
excited state (detector 2) considering (a) 10, (b) 13, and (c) 16 field modes. Both detectors are sharply switched on at the 
same time, and the distance between them is such that the light- crossing time from one to the other is rc- The vertical lines 
represent the time at which this mutual influence between detectors should becomes possible. In each plot we consider what 
happens when we include a different number of field modes in the calculation. Detector 1 begins its evolution in its ground 
state. The solid (blue) lines represent the excitation probability of detector 1 when detector 2 is also started in its ground 
state, whereas the dashed (red) lines are when detector 2 is started in an excited (squeezed with r = 5) state. As expected, 
the initial squeezing of detector 2 contributes to the subsequent excitation of detector 1. However this influence should not be 
able to reach detector 1 until r/rc = 1. By considering different numbers of field modes we explicitly observe the emergence of 
causal signaling as the number of field modes is increased. The parameter values are L = 2it, A = 1/100 and f2 = 9/2 (resonant 
with the 9th mode). 



fects of the other via propagating excitation in the field. 
We will choose sharp switching functions for both of the 
detectors. We will furthermore take the initial state of 
one of the detectors, say detector 2, to be highly excited, 
specifically a single mode squeezed state with a covari- 
ance matrix of the form 



0-2 



,2r 







-2r 



(68) 



where r is the squeezing parameter. Aside from this we 
take the other detector and the field to be in their vacuum 
states. The entire initial detectors+field covariance ma- 
trix is therefore equal to the identity except for the two 
diagonal entries corresponding to detector 2. We then 
switch on both detectors at time r = and compute the 
excitation probability of detector 1 via Eq. (67). 



We place the two detectors at positions xi = L/A and 
X2 — 3L/4, such that they are a distance L/2 apart. 
Since we are taking c = 1 this means that the time from 
the initial switching required for the detectors to become 
in causal contact is Tc = L/2. In Fig. |3]we display the 
excitation probability of detector 1 as a function of t/tc, 
where we show the results including 10, 13 and 16 field 
modes. 

The vertical line in Fig. [3] represents the moment 
in time r = Tc at which the two detectors come into 
causal contact. In each of the plots the solid blue curve 
is the detector 1 excitation probability for when detec- 
tor 2 is initially in its ground state r = 0, whereas for the 
dashed red curve we initialize detector 2 in a squeezed 
state with squeezing parameter r = 5. As expected, 
we observe increased excitations in detector 1 that are 
caused by the propagating field quanta emitted from the 
squeezed detector 2. We see however that if one does 
not include enough field modes then the additional ex- 
citation occurs before the two detectors should be in 



causal contact, therefore implying super luminal signal- 
ing. It is only when we include enough field modes that 
we find the two curves diverging only at r = r^, imply- 
ing that this is when they have mutual influence. Note 
that, of course, even with 16 modes the signaling can be 
seen to be slightly acausal. Further increasing the num- 
ber of modes further improves the causality, although it 
quickly becomes the case that one must include many 
more modes to see a slight improvement. It is only in 
the limit of infinite modes that we have causality in the 
strict sense. 

As a final word on this, we point out that examin- 
ing this emergence is trivially easy using our method. 
First, the symplectic evolution S(t) can be solved an- 
alytically since we use stationary detectors with sharp 
switching functions. Second, in each individual plot from 
Fig. ([3]) the two different curves, solid and dashed, are 
computed using the same transformation S(t). The only 
difference is the initial state <t(0) that we evolve via 
cr(T) = S(T)(T(0)S(r)-^. Here we have only briefly exam- 
ined the emergence of causal signaling — our formalism is 
very well suited for a more complete and encompassing 
study of the effect. Such a study is of great interest to 
the authors and may appear in further work. 



Non-perturbative Detection of the Unruh Effect 
and Thermalization 



In this section we will use the formalism developed 
above to study, beyond perturbation theory, the Unruh 
effect inside a cavity. For a detector with uniform ac- 
celeration a, the worldline (i(r), ^(t)) to be used in our 
Hamiltonian (|64|) is given by 



t(r) = a ^sinh(aT), x{t) = a ^(cosh(aT) — 1), (69) 
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such that the detector is at position a- = at proper time 
r = [36]. 

The excitation stimulated by insufficiently smooth 
switching functions can be somewhat of a problem in 
attempting to observe the Unruh effect in cavities. This 
is because for a given acceleration the crossing time of 
the detector from one side of the cavity to the other may 
be small enough to induce significant noise and wash out 
the Unruh noise. Increasing the acceleration (and there- 
fore the Unruh noise) does not fix the problem because 
then the integration time is even less and that means an 
increased amount of stimulated fluctuations. One way 
of overcoming this difficulty would be to simply increase 
the length L of the cavity, however doing so also increases 
the number of significant modes that must be included 
in the integration and so greatly increases the computa- 
tional effort required. 

For this reason we have opted for a simpler solution. 
Namely, rather than using reflecting (Dirichlet) boundary 
conditions as would be the case in a linear cavity we will 
instead use periodic boundary conditions, as would be 
the case in a periodic waveguide |14j ■ In this setting it is 
physically acceptable for us to arbitrarily increase the in- 
teraction time to values large enough that the stimulated 
noise becomes negligible and thus allowing a clearly ob- 
servation of the Unruh effect. Although in taking this ap- 
proach we are forced to include more modes (because we 
must now include the zero and negative frequency modes 
of the field), this is still more efficient than increasing the 
length of the cavity. Whatever method we utilize, consid- 
ering long interaction times provides the desired results. 
Note that the problem can still be relatively challenging 
from the computational point of view since for extremely 
large accelerations the coupling matrices w(t) and g(T) 
from Eq. ( 46 1 become highly oscillatory very quickly. 



As an aside, there is further reason to use periodic 
boundary conditions over reflecting when studying ac- 
celerated detectors: In a reflecting cavity an accelerated 
detector will observe the field modes becoming increas- 
ingly blue shifted as it travels faster. Eventually even the 
fundamental mode will be shifted beyond the detector's 
resonance frequency, meaning that it no longer resonates 
with any of the modes, at which point the detector effec- 
tively decouples from the field. If we instead use periodic 
boundary conditions then this effective decoupling does 
not occur. To see why this is so, consider a detector 
accelerating to the right. According to this detector the 
left-moving modes of the field undergo an increasing blue 
shift as in the reflecting cavity. The right-moving modes 
however experience a red-shift, meaning that the detec- 
tor resonates with higher right-moving modes over time 
rather than lower. Thus as long as we include enough 
fleld modes in our calculation the detector will continue 
experiencing resonance throughout its evolution. Using 
periodic boundary conditions we therefore avoid the blue- 
shift induced decoupling. This is clearly better for study- 
ing the Unruh effect inside cavities since it is more similar 
to the physics of free space. 



In order to test the Unruh effect with our model we 
must ask two questions. First, after the interaction is 
complete is our accelerated oscillator in a thermal state? 
Second, if so, does the temperature depend linearly on 
the acceleration? To answer the first question let us re- 
call the definition of a thermal state in the Gaussian for- 
malism [35]. In the case of a single- mode thermal state, 
the covariance matrix is diagonal: cr^^^™^ = diag(i^, i'), 
where v is the state's symplectic eigenvalue. The excita- 
tion probabilities in this case follow a Boltzmann distri- 
bution 



„thcrm 



1 



1 



1 



with corresponding temperature 

T ^Vt/\n{l + 2/{v -1)) . 



(70) 



(71) 



As explained in appendix the time evolution gen- 
erated by Eq. (64 1 is given by a multi-mode squeezing 
unitary S. The means that the detector-t-field state af- 
ter evolution is a multi-mode squeezed state of the form 
510). It is known that the reduced state corresponding 
to a subsystem of a multi-mode squeezed state is not 
in general given by a thermal state. Rather, it is given 
by a squeezed thermal state [TT]. This is because the 
squeezing operation S generally includes both inter-mode 
squeezing (which produces thermal subsystems) as well 
as single-mode squeezing. 

The following question then arises: in our specific sce- 
nario does the detector evolve into purely a thermal state, 
or has it also been squeezed? To answer this we compute 
the symplectic eigenvalue v of our detector state cr^, as 
given by the absolute value of either of the eigenvalues of 
the matrix ificr^, and compare the probability spectrum 
of this state with that of a thermal state that is given 
the same symplectic eigenvalue [25] . From the cases that 
we have examined it appears the answer is that while 
the detector does not become exactly thermal, it is very 
nearly so. That is, the multimode squeezing between 
the detector and the field modes is much greater than 
the single-mode squeeze undergone by the detector. To 
conclude this we c omputed the probability of no excita- 
tion from Eq. ( 67 1 as well as p^^'^'^™ and from Eq. 
(70), i.e. what the probabilities of zero and first excita- 



tion would be were our state a thermal state with the 
same symplectic eigenvalue. For small temperatures and 
therefore small excitation probability, which is what we 
will consider here, a good test for thermality is whether 
or not we have po — p^^°™ <C p^^™'^. If this is satisfied 
then the detector is very nearly thermal. For the cases we 
have examined we have found that po ~ p^^''™ is approx- 
imately seven orders of magnitude less than p^^°™, and 
we therefore conclude that our model does indeed pro- 
duce the detector thermality expected. With this first 
question answered we are free to examine the tempera- 
ture dependence on acceleration, where the temperature 



is given by Eq. ( 71 1 
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FIG. 4. Temperature (after evolution) of an accelerated de- 
tector as a function of its acceleration. We observe a linear de- 
pendence as expected from the Unruh effect. The parameters 
used were L = A-k, A(r) = A exp(-r2/2(5^) with A = 1/100 
and 5 = 8/7, and the detector gap was S7 = 4 (resonant with 
the 8th field modes). 



In the example we will present we start with a negative 
initial value for the proper time, < 0, and evolve to 
the positive time Tf = —Ti. This lets us use a long inter- 
action time while minimizing the computational effort. 
That is, we imagine that the detector is injected at high 
velocity into our waveguide such that the acceleration is 
in the opposite direction to its motion. The detector then 
slows down, comes to a stop at r = and then begins 
looping around in the other direction before exiting at 
the same speed it entered with. Again, we are using pe- 
riodic boundary conditions for the field so that this setup 
makes physical sense. We use a Gaussian switching func- 
tion A(r) = Aexp(— r^/2(5^) with i5 large enough that the 
switching stimulation is negligible. It is after this evolu- 
tion is finished that we compute the temperature of the 
detector. 

For a given set of parameters (see caption) we plot in 
Fig. ([4]) the temperature of the detector as a function of 
its acceleration. We see that the detector temperature in- 
deed depends linearly on its acceleration, in concurrence 
with the prediction of the Unruh effect. This is excellent 
since the question of the response of an accelerated detec- 
tor inside a cavity has previously been unexplored, and 
indeed has generated some debate in discussion. Our 
result represents, to the knowledge of the authors, the 
first confirmation of the Unruh effect occurring inside of 
a cavity, and moreover, in a nonperturbative fashion. 

One may worry that extrapolating our data backwards 
it seems that the temperature does not vanish at a = 0. 
This is due to a couple of factors. First, since we are in 
a cavity we should not expect the Unruh effect to hold 
for very low accelerations. For very low accelerations the 
characteristic length c^/aoi the acceleration will be much 
greater than the length of the cavity, at which point we 
expect to see significant border effects, so the response of 
the detector will not be necessarily thermal. Second, as 



one goes to very low temperature the corresponding prob- 
ability of excitation is exponentially suppressed. This 
means that for very small temperatures the Unruh ef- 
fect will be washed out by the switching noise even if the 
switching function is very smooth. 



D. Vacuum Entanglement Harvesting 

In addition to the Unruh effect, another relativistic 
quantum phenomenon that is of great interest is the ex- 
traction of entanglement from the vacuum field. That 
is, two detectors can become entangled by each inter- 
acting locally with a quantum field, even if they remain 
spacelike separated [TH [^TJ [231 [31] ■ Of course it is well 
known that no local operations can increase entangle- 
ment between two quantum systems I37j . In the case 
at hand however there is already entanglement present 
in the vacuum state of the field between spatially sepa- 
rated degrees of freedom, and so by interacting with the 
field locally, multiple detectors can extract this entangle- 
ment to become entangled themselves. This is true even 
if the detectors remain spacelike separated throughout 
their evolution, meaning that they can become entangled 
without any direct causal mutual influence. 

In our detector model we can consider multiple detec- 
tors very easily by simply adding their respective field 
interactions into the coupling matrices w and g used in 
Eq. (46). Once the evolution has been solved and the 



detectors-|-field covariance matrix cr obtained, the multi- 
detector covariance matrix (Td is obtained by deleting the 
rows and columns corresponding to the field. In the case 
of two detectors we obtain a 4 x 4 covariance matrix that 
can be arranged in the form 



f 1 C12 

T 



(72) 



where cri and (T2 arc the 2x2 covariance matrices of the 
detector 1 and detector 2 subsystems, and are of the same 
form as Eq. (65 1. (T12 is a 2 x 2 matrix that provides the 
correlations between the two detectors. For a two-mode 
Gaussian state such as this the entanglement between the 
detectors, as measured by the logarithmic negativity, is 
given by [3S] 



En = max(0, — logv^), 



(73) 



where 



(A ± a/A^ — 4det(Trf) /2 are the symplectic 
eigenvalues of the partially-transposed covariance matrix 
and A = detcri + det(T2 — 2dct(Ti2. 

Using this, we plot in Fig. ^ an example of the loga- 
rithmic negativity between two detectors as a function of 
t/tc, where tc is the time it takes for the detectors to 
come into causal contact. We used sharp switching func- 
tions for both detectors, as we observed that the stimula- 
tion incurred by switching aided in generating entangle- 
ment. Note that due to this sharp switching we needed 
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FIG. 5. Extraction of vacuum entanglement with two detec- 
tors. Notice that the entanglement harvesting starts close to 
the time t/tc — 1, where the detectors come into causal con- 
tact, but still before they reach this time. The parameters 
used are L = 2tt, X = 1/100 and Q = 9 (resonant with the 
18th mode). 



to include many modes, up to 100, in order to find con- 
vergence of our solution. The positions of the detectors 
in this case was chosen to be = L/A and X2 — 3i/4, 
implying that tc = L/2. We indeed observe that, given 
enough time to locally interact with the field, the two 
detectors become entangled. Furthermore we see in this 
example that entanglement is produced before the detec- 
tors come into causal contact, although just barely. 

As a final note, if we compare a harmonic oscillator 
particle detector with the standard qubit detector, one 
can easily show that to second order in perturbation the- 
ory (leading order in this phenomenon) the only differ- 
ence between the evolution of a single qubit and that of a 
single oscillator is that the oscillator develops off-diagonal 
coherence terms in its density matrix. In view of this fact 
we can hypothesize that oscillators are less efficient at ex- 
tracting spacelike entanglement than are qubits, so the 
use of two level systems may be more appropriate in what 
regards field entanglement harvesting. 



V. CONCLUSIONS 

By applying the Gaussian formalism we have addressed 
the problem of time evolution of a particle detector un- 
dergoing relativistic movement inside of a cavity. With 
this, we are able to tackle arbitrarily multimode time- 
dependent problems and solve them nonperturbatively. 
This is markedly different from the standard Unruh- 
DcWitt model that can generally be only solved per- 
turbatively. Remarkably, the only fundamental change 
between the standard approach and our work is that we 
use a harmonic oscillator to describe a detector rather 
than a qubit. 

In addition to being nonperturbative, the methods we 
have presented lead to a computationally very efficient 
way of solving a great range of problems involving an 
arbitrary number particle detectors coupled to quantum 



fields inside a cavity such that the detectors: 1) can be 
undergoing arbitrary relativistic motion, 2) can have ar- 
bitrary quadratic interaction with the field, 3) can start 
in arbitrary Gaussian initial states, and 4) in the field we 
can include any number of modes and impose arbitrary 
time-dependent boundary conditions. The vast range of 
scenarios that this can encompass are all solved by the 
same number-valued, linear, first-order ordinary differen- 
tial equation. We have the additional advantage that for 
a given evolution we do not need to re-solve the equation 
again if we decide to change the initial state. 

To demonstrate this wide applicability we have ana- 
lyzed a range of different problems of interest in gen- 
eral relativistic quantum field theory, obtaining several 
results. 

One of our most important findings is that an acceler- 
ated harmonic oscillator detector in a cavity experiences 
a thermal response with temperature proportional to its 
acceleration. Namely, we have demonstrated that the 
Unruh effect occurs inside of cavities, a scenario that we 
believe has not been previously explored. We empha- 
size that we obtain the thermal response from the move- 
ment of the detector, as opposed to assuming a different 
quantization scheme for accelerated observers and going 
through Bogoliubov transformations. We thereby show 
evidence pointing towards the universality of the Unruh 
effect: a thermal response proportional to the detector's 
acceleration appears even 1) considering a different model 
instead of a standard UdW detector, 2) inside a cavity 
and 3) non-perturbatively. Also we are able to analyze 
the border effects appearing due to the finite size of the 
cavity. 

We have been able to quantify the effects of sudden 
switching for inertial detectors, determining the excita- 
tion probability in different scenarios. The results ob- 
tained are as expected: greater stimulation is generated 
the sharper the switching function is. 

Remarkably, we also quantitatively addressed the 
problem of recovering causality in the context of fully 
relativistic cavity field theory in which there is a UV- 
cutoff. We have shown that whereas a single-mode model 
strongly violates causality, causal signaling emerges as 
more modes in the field are included. Although rigor- 
ously speaking an infinite number of discrete modes are 
required to recover causality, we can determine how many 
modes are necessary to include in the problem to have 
causal signaling for a given time precision. 

We also analyzed the problem of vacuum entanglement 
harvesting. We have shown that indeed harmonic oscil- 
lator detectors inside a cavity can extract entanglement 
from the vacuum field and can achieve this entanglement 
while they are spacelike separated. 

We close emphasizing again the immense applicabil- 
ity of our model. There are many problems of interest 
that can now be easily addressed, and we hope that our 
work will be taken advantage of in the future to explore 
the problems of relativistic quantum physics. For exam- 
ple this work is perfectly suited to studying the model 
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dependence/independence of the Unruh effect: with our 
formaUsm it is trivially easy from the technical point of 
view to modify the detector-field interaction, as well as 
considering problems involving arbitrary trajectories in 
flat space times as well. Applications to cavity settings 
in curved space times can be also thought of. For in- 
stance it could be easily shown how to translate some 
of the results to cavities close to the event horizon of a 
stationary black hole [31], and it would be undoubtedly 
interesting to analyze what would be the behaviour of 
these quantum systems when they experience a dynami- 
cal gravitational collapse [321 SD] ■ 

Note added: During the final stages of this project, 
we became aware of another research group working 
along similar lines [31]. We both agreed to release our 
completed work simultaneously. 
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Appendix A: Detailed derivation of the calculation 
of the Hilbert space evolution 

We will provide here the full detail of the derivation of 



displacement operators are respectively defined as [57] 



the results presented in section IIIB 



As said in the main text, for a quadratic Hamiltonian, 
time evolution can be expressed in terms of displace- 
ments, squeezing, and rotation unitary operations. In 
particular, the unitary evolution we are looking to solve 
for can then be put into the form 

U{t) = e'^M5(z(r))^(/3(T))i?(<A(T)) , (Al) 

where 7 is number valued and the squeezing, rotation and 



S(z) 



R{(j)) = e'(a') <^a^ 



D{(3) 



(A2) 
(A3) 
(A4) 



where z and (p ^-re matrices, and /3 is a column vector. 
The notation z^ is used to represent the Hermitian con- 
jugate of z, the elements of which are number valued. 
Note that we without loss of generality we can consider 
z to be symmetric because it is only the symmetric part 
that contributes to S{z). Also note that <{> must be a 
Hermitian matrix to ensure unitarity of R{4>). 

The exact form of these transformations can be ob- 
tained nonperturbatively by employing a technique in- 
troduced by Heffner and Louisell [28]. We will utilize 
the polar decomposition of z into a product of a hermi- 
tian and a unitary matrix, which can always be achieved. 
This takes the form 



z = re*« = . 



(A5) 



where r and 9 are hermitian matrices, and the sec- 
ond equality results from the assumed symmetry of z. 
From here we wish to evaluate how such operators evolve 
the ladder operators of our system so that we can de- 
termine their corresponding symplectic transformations 
on the phase space and therefore the covariance ma- 
trix ascribed to, for example, a multimode squeezed 
state. Using the BCH's Hadamard lemma, e^Be~"^ = 
B + [A, B] + [A, [A, B]]/2\ + . . . , it is straightforward to 
obtain 

S'^ (z)aS'(z) = cosh(r)a + sinh(r)e'^a^ , (A6) 
R\4))kR{ct))^e''l'sL, (A7) 
^^(/3)aD(/3) =a + /3, (A8) 

and similarly 

S'^(z)a^S'(z) = cosh(r'^)a^ + sinh(r'r)e-'^^a, (A9) 

R\c{))a< RicP) = e-'"t'''k\ (AlO) 

i3^(/3)a^i)(/3)=a^+/3*, (All) 

For our purposes throughout the rest of this appendix 
we will ignore the contribution from the displacement 
operator. This is because we typically interested here 
in Hamiltonians that are quadratic and without linear 
terms. Since linear terms are what drive displacements 
we need not consider them here. Generalizing to include 
displacements is however quite straightforward. 

We will consider an interaction Hamiltonian in the in- 
teraction picture that takes the form 

= (a^)Tw(T)a + (a^)Tg(r)a^ + aTg(r)«a (A12) 

and we wish to solve for the unitary evolution tj (a^, a, t) 
that it generates. 
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Consider now this unitary operator in its normal or- 
dered form [/'"•'(a^, a, r) = [/(a^,a, r), where for exam- 
ple (ada^)^"-* — ^d^d + 1- As explained in (2^ we can 
equally well represent this using a number-valued func- 
tion corresponding to U^"^^ of the form {/("^(a^, a, t) — 
where a. and a* are taken to be col- 
umn vectors consisting of real, independent variables. 
That is, we put IJ into normal form and replace a and 
a^ by vectors of number- valued entries a. and a* . In 
this representation Schrodinger's equation idrUij^ = 
HjIsl' ,aL,T)U{T) becomes 



OT 



where -ffj"^ {a*,a + d/ da* , r) is obtained by putting Hj 
into normal ordered form (which in our case it already 
is) and replacing a and a^ by a -I- 8/ da.* and a* respec- 
tively. What we now have is a set of coupled, ordinary 
differential equations. An ansatz for the solution that 
we will use is C7^"^ = e'^^"*'"''^-', turning the equation 
into one for G. Once the solution has been found we 
can then obtain the normal ordered unitary by replacing 
back a and a* by a and a^ and applying the normal or- 
dering operator C/'"^ = ■.e'^^^ ,a,T). where, for example. 

Following the prescription of [ _2B1 we now want to find 
the evolution eq uatio n of the number-valued function 
[/("). From Eq. we have 



in) 



+ ^^:T]m-\a*a,T), 



da* 



(A13) 



[a* 



w o; 



da 



a 



g 



a 



d 

a(a*)T 

and making the ansatz [7'^"^ = e 
for G: 

.dG , . dG 



+ ^)]^^"^ (AM) 
^ we have the equation 



* 57 = Kfwa + {a*f^-^^ 4- {a*f^a* 



a^ ^a 



^T„H 

a g 



dG 



dG 



da* d{a* 



dG 



dG 



d 



Tg 



dG 



d[a*y^ da* d{a*y^ da*' 
Additionally we can make the educated ansatz 
G = {a*fTi(T)a + (a*)'^C(r)a* + Q;'^F(r)a 



(A15) 



-a(T), 
(A16) 



where £), C and F are matrices. In general we should 
also include terms linear in a and a* corresponding to 
phase space displacements, but in our case they will be 
absent due to the lack of linear terms in the relevant 
Hamiltonians and so we will not consider them. From 
here it is easy to show that 

dG 

— =Da + 2C,a*, (A17) 



where = (C -I- C'^)/2 is the symmetric part of C. 
The transposed version of this relation follows trivially. 
Lastly, it is easily shown that 



d 



d{a*Y 



g«|^=2Tr(gHa 
da* 



(A18) 



Given these relations it is now a simple matter of com- 
paring coefficients between the right and left sides of Eq. 
(A15I. Doing so, we find the coupled set of differential 



equations 



lA = 2Tr(gHa), 

it = 4C,g«C, + 2wC, 

iti = (4C,gH+w)(D4 



g, 



I)g"(D- 



(A19) 
(A20) 
(A21) 
(A22) 



where I is the identity matrix and we have initial condi- 
tions ^(0) = and C(0) = D(0) = F(0) = 0. 

These equations can be numerically solved, although 
we will find that for our purposes the only one that actu- 
ally needs to be solved is the equation for C. This is be- 
cause C fully determines the squeezing matrix z — re'^, 
which since our system is initially in the vacuum state 
is all that we need (since the vacuum is invariant under 
rotations). That is why we can thus ignore the rotation 
and effectively set cp — 0. For a more general initial state 
one would need to additionally solve for D in order to 
compute 0. Note also that one will never have to solve 
Eq. (A22| for F; it can be expressed purely in terms of 
C and D and is therefore a redundant variable. 

Note that the form of these equations are entirely in- 
dependent of the specific coupling matrices w and g that 
we choose. We are therefore free to choose an entirely 
different interaction Hamiltonian and the evolution will 
still be represented by these equations. Once solutions 
have been found we can return the (normal ordered) uni- 
tary to its operator form via U^"^ = -.e'^^'^ , which 
from Eq. (A16) gives us 



^F(r) 



(A23) 



Our problem is thus essentially solved; the last task 
required is to overcome the normal ordering and to put 
this unitary into a form that we are familiar with. The 
procedure for doing this is given in |27j , but we will reit- 
erate it here in somewhat more detail. We know that U 
should be a product of rotation and squeezing operators, 
along with possible phases: 



U = e*''S'(z)i?(0), 



(A24) 



where S and R are given by Eqs. (A2 A3). Note that 



in our case U contains no displacement operator because 
Hi contains no linear terms. We know that the solutions 
Eqs. (A23) and (A24| must be equivalent, and the task 
is now to find the relation between 7, z, and (p ^-^id the 
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results obtained for A, C, D and F. Trivially we see along with the fact that F is symmetric, as can be seen 



that ij = A and from Eq. (A19), however this merely from Eq. (A22). With these, Eq. (A27) gives us 



contributes an overall phase to the evolution and we will 
therefore ignore this contribution henceforth since it does 
not contribute to the physics 

In order to find z = re^^ — e*^ r""" and <p we recall the 
action that Eq. (A24| will have on the ladder operator: 



U'^aU = cosh(r)e^'^a + sinh(r)e*^e~'<^^a^ (A25) 
U^a'U = cosh(r^)e-'<^^a^ + sinh(r'^)e"*^^e"^a. (A26) 

We now need to compute tj^ktj and tj^s^tj from the uni- 
tary in Eq. ( A23 1 in order to compare. To this end we use 
the identities [a, tj] = d(j /ds^ and [a^, tj] = ~dU /da, or 
equivalently 



U'faU = 



dU 



i, U'^a^U = -U^ 



dU 



it 



(A27) 



da^ ' ^ 9a 

To evaluate the right-hand sides we use the identities: 
d 



„(at)TDa. ^ :e(«^)^i3«:D£ 



9a 

d 

(/ + DT)a^e('^')"°=^ 

^^tga^Fd ^ ga^Fa^^^t _ 2Fa), 



e('^')"°-at, 



(A28) 
(A29) 
(A30) 



U^'aU = [(D + I) - 4C,(D'^ -t- l)-^F]a 

+ 2C,(D^ + I)-la^ (A32) 

U^a^U = {T>^ +I)-\-2Fa + a'). (A33) 

Since these equations are just the adjoints of each other 
we are able to compare the two and determine the addi- 
tional relations 



F = -(D* +I)-iC*(D + I) 
I-4C,C: = (D + I)(D^+I) 



(A34) 
(A35) 



Given all of this, we find indeed that (/''aU and WaW 
are of the form given in Eqs. ( A25|A26 1 where we identify 



C, = -tanh(r)e^^ 
D + I = scch(r)e*'^. 



(A36) 
(A37) 



Thus, once we have integrated Eqs. ( A20|A21 ) for C and 
D we can use this result to solve for the corresponding 
squeezing and rotation matrices r, 9, and (p and, via Eq. 
(61 1, obtain the covariance matrix in which all properties 



(A31) 

of the final state are encoded. 
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